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Abstract 

D ■ We review the present knowledge on the heavy quark interaction. The 

1^ . framework is the NRQCD effective theory and the interaction is expressed 

VO \ in terms of Wilson loop chromoelectric and chromomagnetic insertions. 

> 1 Introduction 

o 






X 



The experimental evidence that for heavy quark bound states like 66, cc, ... all 



CN \ the splittings are considerably less than the masses suggests that all the dynam- 

qT \ ical energy scales of these systems are small with respect to the quark masses. 

0> ' As a consequence the quark velocities v are small and these systems can be con- 

^ ' sidered as non relativistic. The hierarchy of the scales is then the typical one 

'^i of a non relativistic system. Called m the mass of the heavy quark, the quark 

Qh! momenta scale like mv and the quark energies like mv'^. Trickier is the situation 

(-H \ in the gluon sector, but the binding interaction is essentially characterized by 

the same energy scale distribution. Therefore the dominant gluon interaction 

among heavy quarks appears "instantaneous" . A potential picture should hold, 

at least in first approximation, and the energy levels can be obtained by solving 

5^ \ the corresponding Schrodinger equation. In particular for infinitely heavy quarks 

the spin splittings vanish and the system can be described only by means of a 

single static central potential. 

Indeed, potential models have been quite successful [jl| in explaining the 
hadron spectrum. However, being these phenomenological models, their con- 
nection with the QCD parameters is hidden, the scale at which they are defined 
is not clear and they cannot be systematically improved. In this paper we outline 
the kind of rigorous description of quarkonium that can be obtained from QCD 
at the present. 

Several perturbative evaluations of the quarkonium potential have been per- 
formed in the last twenty years and recently the complete ttg corrections to 
the levels have been calculated [@. The main difficulty in these calculations is 
connected with the inclusion of nonperturbative contributions which not only set 
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a bound on the reachable precision, but turn out to be relevant for quantitative 
predictions even on the bottomonium spectrum. A typical approach is to consider 
all the nonperturbative physics to be encoded into few local condensates. This 
corresponds to assume that the typical length associated, for instance, with the 
nonperturbative gluodynamics is larger than any other length scale of the system 
and can be put equal to infinity. In a pure analytic calculation this statement 
is also dictated by the necessity to have a small number of nonperturbative pa- 
rameters and maintain predictability. Even for the ground state this assumption 
is doubtful, but surely it does not hold for large radii quarkonia (i.e. excited 
heavy mesons) where the nonperturbative gluonic length cannot be considered 
large with respect to the size of the bound state. Therefore a proper treatment of 
nonperturbative effects includes non-local condensates. Another way to say the 
same thing is to notice that the contribution to the levels associated with the lo- 
cal gluon condensate is proportional to n^ {agF"^ (0)) . By increasing the principal 
quantum number n beyond the ground state, it grows very soon out of control. 

A way of including nonperturbative effects, as they are, in the evaluation of the 
quarkonium spectrum is lattice QCD. This technique is becoming more and more 
successful and in the near future it is expected to be the only competitive one (see 
for instance [Q and references therein). There are different ways in which lattice 
QCD calculations can be performed. We mention lattice NRQCD f^ where the 
quarkonium spectrum is directly evaluated on the lattice by means of an effective 
action derived from QCD by an expansion in the quark velocities. Very close is the 
approach we will discuss in the following where the QCD Lagrangian is replaced 
by its effective non relativistic formulation before doing any lattice evaluation. 
All the non perturbative physics is encoded in this way in the Wilson loop made 
up by the quark trajectories and in field strength insertions on it. These (non 
local) objects are then evaluated on the lattice. We call this the Wilson loop 
approach. The advantage is that in this way our expressions are safer to handle 
for lattice purposes, since less affected by finite size effects, and easier to treat also 
for analytic purposes like the implementation of vacuum models. In particular 
we get an expression for the heavy quark-antiquark potential. In the limit where 
the insertions of two field strengths on the Wilson loop can be approximated by a 
local condensate, one gets back the "improved" perturbative expression discussed 
previously. 

The (Wegner)-Wilson loop formalism has a long story. It was first suggested 
by K. Wilson that the object called after him Wilson loop would be the relevant 
one in order to describe confinement. The strong coupling expansion suggested 
an area law behaviour further confirmed by lattice simulations, which were born 
with this pioneering work. In [^ it was shown how to relate rigorously the static 
Wilson loop with the static quark-antiquark potential. On the same line few 
years later also spin dependent corrections to the potentials where expressed in 
terms of Wilson loop and chromoelectric and chromomagnetic field insertions ||^. 
In particular it was proven that a nonperturbative behaviour in the static poten- 



tial must give rise to nonperturbative spin dependent corrections as well. Non 
spin dependent corrections were treated in the same framework some time later 
0. This was the situation pictured at the beginning of the '90 in [Q]. Several 
problems were still open. No relativistic formulation in terms of the Wilson loop 
was available. There was an apparent mismatch between the Eichten-Feinberg- 
Gromes expression for the spin dependent sector of the potential (analytic in 
the quark masses) and the perturbative one loop estimate (containing logarithms 
of the quark masses). The inclusion of non potential terms was not understood. 
While not all these problems have been solved, some remarkable progress has been 
achieved in the last years. Attempts in the direction of a relativistic formulation 
have been done in |jlO|]. In [|ll] it was shown that, by performing properly the 



matching between the effective theory and QCD, Wilson coefficients carrying log- 
arithms of the quark masses appear in the Eichten-Feinberg-Gromes expression 
for the potential. The inclusion of non potential contributions in this frame- 



work is now in progress |T^ and the way to precision calculations in quarkonium 
seems finally to be open. Here we only mention that this goal is met by a new 
effective theory built from NRQCD, where explicitly potential and non potential 
terms have been separated [0. Finally, all the Wilson loop averages relevant 
for the potential have been calculated on the lattice and the bottomonium and 
charmonium spectra have been calculated with good agreement with the data 

ini,[T5ii. 



In the meantime there have been several progress in the building of a non 
relativistic effective theory from QCD, mainly due to the success of Heavy Quark 
Effective Theory in describing heavy-light systems |T^. Wilson coefficients have 



been calculated at higher order and the role played by reparameterization invari- 
ance has been better understood. 

Here, we summarize the present level of understanding of the heavy quark 
potential which contributes to the energy levels of quarkonium at the order v^. 
The framework is NRQCD, the tool the Wilson approach. 

2 NRQCD 

In order to define an effective theory we typically need three ingredients: an 
effective Lagrangian, a regularization scheme and therefore a matching scale and 
a power counting set of rules. The effective theory we will use is NRQCD ||^. 
Let us discuss its key ingredients. 

The NRQCD Lagrangian is obtained from QCD by expanding with respect 
to the heavy quark masses. The matching with QCD is performed like in Heavy 
Quark Effective Theory (HQET) [^, |18|. We emphasize that in order to build 
up the effective Lagrangian from the QCD Lagrangian, we have to ignore the 
specific dynamical problem we are dealing with and expand with respect to the 
heavy quark masses which are explicit parameters of QCD. Typically the effective 



Lagrangian turns out to be the sum of a pure gauge part Lg plus two, four, ... 
fermion terms (L2f, -L4f, ...). Since we are interested in two-body bound states 
we will take into account only two and four fermion terms (terms involving more 
fermions will contribute only in intermediate states). Therefore the NRQCD 
Lagrangian obtained from QCD by expanding with respect of the mass mi of a 
heavy quark and the mass m2 of a heavy antiquark is given by [1T7| , |18 
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Qj are the heavy quark fields and the covariant derivative is defined as D^ = 
d^ + igA^jT"-. The explicit form of L4f is given in |T^. Lg is the SU(3) Yang-Mills 
Lagrangian modified in order to give rise to an effective ag running with 2 (heavy) 
flavors less. 

The effective Lagrangian (|I|) is not renormalizable. Therefore, it is necessary 
to regularize it. In a given regularization scheme the reproduction of the correct 
ultraviolet regime of QCD is obtained by means of the Wilson coefficients. The 
effective Lagrangian is complete only once these coefficients are given. The Wilson 
coefficients are evaluated at a matching scale where perturbation theory still 
holds. They encode the ultraviolet regime of QCD up to a given scale /x order by 
order in the coupling constant Og. Renormalization group (RG) transformations 
should be used in order to resum leading log contributions (~ In m/fi). The 
matching coefficients are known in the literature at different level of precision (i.e. 
at different order in the coupling constant), see [p!7| , |1^, |T^ Here we remember 
only that reparameterization invariance |T^ (i.e. the invariance of the effective 
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Lagrangian with respect a variation of v) fixes q" 

Heavy quark bound states are characterized by a dynamical dimensionless pa- 
rameter, the quark velocity v. As explained in the introduction, this parameter is 
small. From phenomenological potential models [Q and from lattice studies |]1^ 
we get the usually accepted values of (f^) ~ 0.07 for the bottomonium system 
and (f^) ~ 0.24 for the charmonium system. This allows the ordering of the 
energy scales of the problem, m, mv and mv'^. The first scale has been explicitly 
integrated out in the Lagrangian (ffl) and, as discussed above, the contributions 



coming from it are carried order by order in ctg by the Wilson coefficients. The 
last two scales, sometimes known with the name "soft" and "ultrasoft" respec- 
tively, are still mixed up. The relation between v and the QCD parameters is 
in general unknown. It is expected to be the result of perturbative and nonper- 
turbative effects. For infinitely heavy quarks v coincides with as like in QED 
where, for instance in the hydrogen atom, the hyperfine structure constant a is 
equal to the electron velocity. The main point is that the existence of this small 
dynamical parameter enables us to set up power counting rules for the opera- 
tors in the effective Lagrangian. For the sake of simplicity we reproduce here 
the rough power counting argument of 0], keeping in mind that an exact power 
counting cannot be done until soft and ultrasoft degrees of freedom have not been 
completely disentangled. Noticing that the number operator for heavy quarks, 

/ d?xQ'^{x)Q{x)^ does not depend on f, one gets Q ~ (?7if)^/^. The kinetic en- 

ergy / d?xQ\x) Q{x) is by definition of order mv"^ and therefore D ~ mv (i.e. 

the contribution of D to the levels is of order mv) . Finally, using the equation 
of motion one gets the other counting rules Dq ~ mv"^, qAq ~ mv'^ gA ~ mv^, 
gE ~ rn^v^ and gB ~ m?v'^. With respect to these rules the Lagrangian of Eq. 
(H) is accurate up to order 0{v^) of the levels (or up to 0(f ^) with respect to 
the leading contribution). This should guarantee a rough 10% of accuracy on the 
absolute value of the levels. 

Concluding this section, we stress that the NRQCD power counting defined 
above is not the same as that one used in HQET. In particular, in HQET the 
kinetic energy D^/2m is suppressed by a factor Aqcd/''^ with respect to the 
operator Dq while in NRQCD the two operators are of the same order. As 
a consequence the heavy quark propagator contains in NRQCD a kinetic part 
which is absent in HQET where the heavy quark propagator is static. In a very 
general way one can say that these differences are due to the fact that, even if 
the effective Lagrangian is essentially the same, the physical problem is different 
and it is the physical problem which defines the counting rules. 

3 The Wilson loop formalism 

The next step is to derive the heavy quark interaction in the so-called Wilson 
loop formalism. In this context the use of the effective Lagrangian (Q) (with tree 
level matching) was first suggested by L. S. Brown and W. I. Weisberger, later 
by E. Eichten and F. Feinberg in [0, §]. A one loop RG improved calculation 



was done in |]rT|. In those papers an expansion is performed around the static 
solution. Here, we adopt the approach of ||^ where the kinetic energy was kept 
during all the calculations. 

The 4-point gauge invariant Green function G associated with the Lagrangian 



(|1|) is defined as 

G(Xi,|/i,X2,l/2) = (O|(55(x2)0(x2,Xi)Qi(Xi)(5l (1/1)0(1/1, |/2)Q2(y2)|0), 

where 0(x2, xi) = exp < —ig / ds (x2 — xi)''A^(a;i + 5(0:2 — a;i)) > is a Schwinger 

line added to select the singlet state contribution. For large time separations the 
string vanishes. After integrating out the heavy quark fields Qj and Qj, G can 
be expressed as a quantum-mechanical path integral over the quark trajectories 
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where Zj = Zj{t) and 1/2 = Z/i — ~T/2, X2 = x^ = T/2. The function U describes 

the heavy quark interaction. Since the kinetic energy has been properly isolated, 

it is consistent with the counting rules given in the previous section to expand the 

interaction around the static quark trajectories Zi = (t, r) and Z2 = {t, 0). More- 

1 fT/2 
over, we define the heavy quark-antiquark potential as \^ = lim — / dtU, 

T^oo T J -T/2 

with the warning that, having soft and ultrasoft degrees of freedom not been dis- 
entangled in NRQCD, non potential terms could still contribute to some extent 
ioV. 

Working out this expression |^, |T^ we get 

V = Vo + spin indipendent terms idss + Cp d^s) (5^(r) 

mim2 
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S(-5') = <t'^-'^/2 and L'^-'^ = r x pj are the spin and orbital angular momentum 
operators respectively. The matching coefficients dg come from the 4-fermion 
sector L^i |TB| and 20^^. = Cp'g ± Cp\. The static potential Vq is 



V^{r) = hm ;^ln(iy(ro)), l^(r) = P exp i-tg i dz^A,{z)\ . (4) 



where F is the loop made up by the quark trajectories zi and Z2 and the endpoints 
Schwinger strings, and the static loop To is a r x T rectangle. The bracket 
means the color trace of the average over the gauge fields weighted by the gluon 
Lagrangian Lg. The "potentials" Vi, V2, ... are scale dependent gauge field 
averages of electric and magnetic field strength insertions in the static Wilson 
loop, see [0, |1^ . 



Since the spin independent corrections come from the terms Dq and D^/2m 
of the Lagrangian (Q) whose matching coefficients are protected by reparame- 
terization invariance, they are scale independent. An evaluation of the so-called 
momentum dependent corrections can be found in ||^. Terms involving logarithms 
of the quark masses are present. As a consequence of the matching procedure 
they are all encoded in the matching coefficients. Therefore, as argued in [ p!l[| , 
the correct handling of the matching allows finally the agreement between the 
potential derived here evaluated in the perturbative regime with the traditional 
QCD one loop perturbative calculation [|, |21|. For a careful discussion of this 
point see |T^ . 



The "potentials" are known exactly (up to a given order in ctg) only in the 
perturbative regime, i.e. in the short range behaviour. Nevertheless there exist 
some exact relations between them which hold at any range. For instance, from 
Lorentz invariance it follows that |^ 

Kir) + Vlir) - V^ir) = 0. (5) 

Since for reparameterization invariance c^ = 2cp — 1 |T^, equation (^) holds at 
any scale /i. Similar relations exist for the momentum dependent "potentials" of 

Ref. [01. Moreover the scale independence of the potential (BF) (i.e. /x— — = 0) 

establishes several relations between quantities at different renormalization scales 



0- 

The potentials Vq, Vi, ... are suitable for direct lattice computation and for 

analytic evaluation inside a QCD vacuum model. The first possibility relies on 

the fact that all the dynamical quantities are expressed in terms of field strength 

insertions on a static Wilson loop which is an object traditionally measured on 

the lattice. Such an analysis has been performed in |T^. All the "potentials" 



have been measured on the lattice and the previous mentioned exact relations 
have been used to check the accuracy of the results. Using the parameterized 
form obtained in this way for the "potentials" the bound state equation has been 
solved and the bottomonium and charmonium spectra evaluated. The agreement 
with the experimental data is found to be quite good. On the other hand, an 
analytic evaluation of the heavy quark potential inside a QCD vacuum model 
turns out to be very convenient in this framework since by means of functional 
derivatives all the averages of field strength insertions on the Wilson loop can be 
expressed in terms of the average of the non-static Wilson loop alone. This is 
usually a quantity provided by QCD vacuum models. A study of different models 



in this framework has been done in [|r^, ^ where a comparison with the existing 
lattice data has also been provided. Up to now these are not accurate enough to 
really discriminate between different models. More precise lattice measurements 
will be performed in the near future providing in this way strong constraints on 
all the infrared QCD models with predictability on the long-range quarkonium 
interaction. 

Finally we emphasize that until soft and ultrasoft degrees of freedom have 
not been disentangled an exact value in v cannot be assigned to each term of the 
effective Lagrangian and the power counting has to be interpreted at the leading 
order. For instance the 0(1) NRQCD Lagrangian, L = Qi [iDq + d"^ /2mi) Qi + 
Q2 {—iD^ + d'^ /2ni2) Q2, does not contribute to (H) only with the static potential 

Vq. Since the corresponding Wilson loop P exp < —ig j> dz^AQ{z) > is a function 

of the non-static loop F, its expansion produces spin independent terms of order 
0{v) as well. 

4 Comments and outlook 

In the framework of NRQCD we have shown how to get the complete order 0{v^) 
expression of the heavy quark- ant iquark potential in terms of field strength in- 
sertions on a static Wilson loop. This has several advantages. Nonperturbative 
contributions can be evaluated by means either of traditional lattice calculations 
or of different QCD vacuum models. Having worked out the matching procedure, 
we find that the potential of Eq. (H) is consistent in the short range with the 



existing perturbative calculations and with the lattice data [T^. These are sen- 
sitive to one loop and in some cases to the next to leading correction too. We 
stress here that terms proportional to the static potential in Eq. (0) have to be 
protected from the running. For instance, using Eq. (^) and reparameterization 
invariance, the spin-orbit interaction term ~ 2CpVl + CgV^ can be written as 
2CpV2 — Vq where it is apparent that no Wilson coefficient multiplies the static 
potential. Moreover we can draw the following consequence about Vg': either V2 
does not contain any nonperturbative contribution at all (as present lattice data 
seem to suggest) or, if V2 contains some nonperturbative contributions, they have 

to satisfy the RG equation -rCtVo = 0. The same argument suggests for the 

dfi 

Darwin term that AVj^ = — AVo+ perturbative contributions. This gives some 

constraints on the QCD vacuum models which, introducing at some point some 

approximations, lose scale invariance. 

All the corrections to the effective Lagrangian discussed here are relevant 

in order to obtain the quarkonium spectrum with an accuracy of 0{v^). For 

some Wilson coefficients only the leading as contribution needs to be taken into 

account. This is no more the case if we aim to reach a 10% accuracy in the 

quarkonium spin splittings. Being these quantities an order 0{v^) effect, a 10% 



accuracy is achievable only if 0{v^) and 0{asV^) effects are calculated as well. 
Therefore operators of order 0{v^) should be added to the effective Lagrangian 
[0. The inclusion of such operators in Eq. (^ in terms of (two and three) field 
strength insertions on a static Wilson loop, is only a technical problem, but has 
not be done so far. The main reason is that non potential contributions are 
expected to become even more important. Ultrasoft gluonic degrees of freedom, 
characterized by a time scale l/mv"^, exist. In perturbative QCD they are respon- 
sible, for instance, of the Lamb-shift (which is a agV^ effect). In the language of 
the Wilson loop, this means that the interaction U (see Sec. 3) can be affected 
by non potential contributions which have to be subtracted from V in order to 
define properly a heavy quark potential. Since ctg is not a small parameter at 
the ultrasoft scale, these non potential terms could be an effect of order 0{v^). 
Therefore, precision studies of quarkonium have to take into account it. An ap- 
proach was recently proposed in [^, |T^. The ultrasoft degrees of freedom are 
integrated out directly from the NRQCD Lagrangian giving rise to another ef- 
fective theory, called potential NRQCD, where all the energy scales of the bound 
state are disentangled explicitly. The advantages are enormous, since in the new 
theory potential and non potential contributions are clearly separated. The novel 
feature is that the matching this time takes place in a energy region dominated 
by the nonperturbative physics. 
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